Signatures of homoclinic motion in quantum chaos 



in 
o 
o 

(N 
oo 



Q 
U 



(N 
> 

(N 
O 
(N 
O 

o 



> 
X 



D. A. Wisniacki/'^ E. Vergini,!'^ R. M. Benito,'^ and F. Borondo^'E 

' Departamento de Qmmica C-IX, Universidad Autonoma de Madrid, Cantoblanco, 28049-Madrid, Spain. 
^ Departamento de Fisica "J. J. Giambiagi", FCEN, UBA, 14-28 Buenos Aires, Argentina. 

14.29 Buenos Aires, Argentina. 
28O4O Madrid, Spam. 

(Dated: February 8, 2008) 



J. 

'"^ Departamento de Fisica, Comision Nacional de Energia Atomica. Av. del Libertador 8250, 
Departamento de Fisica, E. T.S.I. Agronomos, Universidad Politecmca de Madrid. 



Homoclinic motion plays a key role in the organization of classical chaos in Hamiltonian systems. 
In this Letter, we show that it also imprints a clear signature in the corresponding quantum spectra. 
By numerically studying the fluctuations of the widths of wavefunctions localized along periodic 
orbits we reveal the existence of an oscillatory behavior, that is explained solely in terms of the 
primary homoclinic motion. Furthermore, our results indicate that it survives the semiclassical 
limit. 
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Periodic orbits (POs) are invariant classical objects 
of outmost relevance for the understanding of chaotic 
dynamical systems. The pioneering work of Poincarc 
demonstrated the importance of these retracing orbits in 
organizing the complexity of the corresponding classical 
motion. Semiclassically, Gutztwiller l] developed his cel- 
ebrated trace formula, able to quantize chaotic systems, 
solely in terms of POs information. Quantum mechani- 
cally, striking manifestations of POs have been reported 
in the literature. One of the most important are "scars" 
an enhanced localization of quantum density along 
unstable short POs in certain individual cigonfunctions. 
This phenomenon was first noticed in the Bunimovich 
stadium billiard and subsequently studied system- 
atically by Heller, who constructed a theory of scarring 
based on wave packet propagation Another impor- 
tant contribution to scar theory is due to Bogomolny , 
who derived an explicit expression for the PO contribu- 
tions to the smoothed quantum probability density over 
small ranges of space and energy (i.e. average over a large 
number of eigenfunctions) . A corresponding theory for 
Wigner functions was developed by Berry Recently, 
there has been a flurry of activity focusing on the influ- 
ence of bifurcations (mixed systems) on scarring . 

The existence of a scar implies a clear regularity in the 
corresponding quantum spectrum, related to the period 
of the PO. In time domain, the dynamics of a packet 
running along a PO induce recurrences in the autocorre- 
lation function, that when Fourier transformed define an 
envelope in the spectrum, giving rise to peaks of width 
proportional to the Lyapunov exponent. A, at energies 
given by a Bohr-Sommerfeld (BS) quantization condi- 
tion XM- 

In this Letter, we demonstrate the existence of an ad- 
ditional superimposed spectral regularity also related to 
scarring. It is originated by the associated homoclinic 
motion and is given by the area enclosed by the stable 
and unstable manifolds up to the first crossing. To un- 
veil this regularity we consider the fluctuations of the 



spectral widths corresponding to localized wavefunctions 
along unstable POs. Our data demonstrate that these 
fluctuations have a surprisingly simple oscillatory behav- 
ior essentially governed by the quantization of only the 
primary homoclinic dynamics. We provide an explana- 
tion of this result in terms of the coherence of this clas- 
sical motion 0, which constitute the natural global ex- 
tension of the local hyperbolic structure around the PO. 
Furthermore, our numerical results indicate that the ob- 
served oscillatory behavior do not vanish as ?i — > 0. 

It should be remarked that, contrary to the previously 
described results on scar theory, our study implies dy- 
namics beyond the Ehrenfest time (^ |ln?i|), when the 
PO manifolds start to cross and the homoclinic tangle de- 
velops, giving rise to subtle quantum interference effects. 
Other interesting papers, also conside ring these longer 
times, are due to Tomsovic and Heller \i_u\, who showed 
how to construct a valid semiclassical approximation to 
wavefunctions past this limit; this came as a surprise 
since the underlying chaos has had time to develop much 
finer structure than a quantum cell. Similarly to what 
happens in energy domain with Gutzwiller trace formula, 
this theory unfortunately requires the use of a number 
of homoclinic excursions increasing exponentially with 
time. Our result indicates that there exists some proper- 
ties of these long term dynamics that can be understood 
in terms of a small number of classical invariants. 

The ideal tool to carry out our investigation are non- 
stationary wave functions highly localized on POs. The 
overlap with the corresponding eigenstates provides in- 
formation on how these structures appear in the spectra, 
and how they are embedded in the quantum mechanics 
of the system. Such functions, associated to a PO, can be 
obtained by using the (dynamical) information up to the 
Ehrenfest time, thus including the hyperbolic structure 
formed by the corresponding unstable and stable mani- 
folds. In this way, the associated wavefunctions spread 
along these dynamically relevant regions of phase space 
(see illustration below in Fig. This can be done in 
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FIG. 1: Configuration space (a), piiase space (b), and wave- 
lengtli spectrum weiglits (c) representations of a scar wave- 
function along tiie iiorizontal periodic orbit of a desym- 
metrized stadium billiard corresponding to BS wavelength 
fees = 211.665. The unstable and stable manifolds of the 
orbit are also plotted in panel (b). In panel (a) the Birkhoff 
coordinates on the boundary of the billiard are shown. 



different ways [lllll^ : in particular, in Ref. Ti'a defini- 
tion of scar functions, based on transversally excited res- 
onances along POs at given BS quantized energies with 
minimum dispersion, is provided. The width of these 
functions in the energy spectrum is given by, 
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where the narrowing factor |ln?i|~^ comes from the in- 
clusion of the hyperbolic structure ^bis respect, 
it should be noticed that although this approach incorpo- 
rates the motion along the manifolds up to the Ehrenfest 
time, interference effects due to the intersections of the 
manifolds are not included. 

The model used in our calculations is the fully chaotic 
system consisting of a particle of mass M confined in 
a desymmetrized Bunimovitch stadium billiard, defined 
by the radius of the circular part, r — 1, and the en- 
closed area, 1 + 7r/4. To calculate the corresponding 
eigenstates, Dirichlet conditions on the stadium bound- 
ary and Neumman conditions on the symmetry axes are 
imposed. We will focus our attention on the dynamics 
influenced by the horizontal PO. The corresponding scar 
wavef unctions, I7), are calculated using the method of 
Ref. [13. 

In Fig. n we show, as an example, the results for the 
scar state with wavelength /cbs = 211.665. This value 
was obtained from the BS rule, k^s = (27r/L/f)(n + 
vh/^), where n — 134 is the excitation number along 
the orbit, Lh = 4 its length, and vh ~ 3 the cor- 
responding Maslov index. In it, all features discussed 
above are observed. Namely, the Husimi based quan- 
tum surface of section jl^ spreads along the manifolds 
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FIG. 2: Relative width for the scar states along the horizontal 
periodic orbit as a function of the Bohr-Sommerfeld quantized 
wavelength (inset), and its Fourier transform. 



structure associated to the horizontal PO, and the state 
appears in the spectrum distributed among the differ- 
ent eigenstates, within a given range around the 
corresponding BS quantized wavelength, with intensity 
weights If, = |(Ai|7)P- 

The energy width of the associated envelope can then 
be defined as (we take M = in such a way that 
is the energy of the particle) 



(2) 



When calculated, this quantity is a decreasing oscillatory 
function of fees , that can be more adequately studied by 
defining its relative variation with respect to the corre- 
sponding semiclassical value. 
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The corresponding results are shown in the inset to Fig.|21 
As can be observed, it exhibit an oscillatory behavior, 
with an amplitude representing only a small fraction of 
(Tsc . When this signal is Fourier analyzed (see main body 
of Fig. 121) it is seen to have only two relevant compo- 
nents, corresponding to the peaks appearing at values 
of the action S = 0.633 and 1.007, respectively. (Notice 
that, in our case, S has units of length since the total lin- 
ear momentum of the particle has been set equal to one). 
The BS quantization of the wavelength implies that the 
Fourier transform of (Troi is a periodic function with pe- 
riod Lh- In this respect, corresponding peaks appear at 
S = -3.367 and -2.993. 

To interpret this result, we consider only the primary 
homoclinic motion corresponding to the horizontal PO. 
This follows the analysis of Ozorio de Almeida j^, who 
studied the quantization of homoclinic orbits, that were 
thought of as defining an invariant torus, approached by 
series of satellite unstable POs. There are two topologi- 
cally distinct types of primary homoclinic points, which 
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FIG. 3: (a) Phase space portrait in BirkhofF coordinates relevant to our calculations. Label A indicates the position of the fixed 
point corresponding to the horizontal periodic orbit. The associated unstable and stable manifolds are represented in full and 
dashed line, respectively. They cross, with different topology, at the primary homoclinic points hi and h2, which map into h'j, 
hi', andh^, h^', .... 

(b) Primary homoclinic areas Shi ^-nd Shj associated to the horizontal orbit. 

(c) Satellite periodic orbits converging to the homoclinic points hi and h2. See text for details. 



appear as the first and second intersections of the asso- 
ciated manifolds. The situation is shown in Fig. |21(a), 
where a picture of the relevant phase space, in BirkhofF 
coordinates, is presented. In it, we have plotted the fixed 
point corresponding to the horizontal PO (labelled A), 
and the emanating unstable (full line) and stable (dashed 
line) manifolds. These manifolds first cross (with differ- 
ent topology) at points hi and h2, and afterwards at the 
sequences h'^, h", . . ., and hJ,, hj, ... (and their reflec- 
tions with respect to p — 0), as fixed point A is ap- 
proached. These infinite sequences of points, which are 
the dynamical images of hi and h2, constitute the pri- 
mary homoclinic orbits, that define two relevant areas in 
phase space, denoted hereafter as S'hi and S'hj [see shaded 
regions in Fig. |2| (b)] , important for the discussions pre- 
sented below. Furthermore, these primary homoclinic 
orbits accumulate infinite sets of satellite POs with in- 
creasingly longer periods, as discussed in Ref. The 
first members of these two families, converging to the pri- 
mary homoclinic orbits are presented in Fig. |31(c). Ac- 
cordingly to our notation, based on the number of times 
that the POs intersect our surface of section, orbits 3~ 
and 3-1- consist of fixed points near hi, h'l (or h2, h^, 
and their reflections with respect to p = 0; orbits 4- and 
4-1- consists of fixed points near hi, h'l, h" (or h2, hj, hj), 
and their refiections; and so on for the longer POs in the 
family. Notice that the refiection with respect to p = 
of points over 9 = or 5 = 1-1- 7r/2 gives the same point. 
Moreover, these POs can be grouped in pairs, in such a 
way that 3- and 3-1- represent the first approximation to 
the primary homoclinic motion (the - indicates the PO 
with shortest length of the pair), 4~ and 44- the second, 



etc. Actually, it is apparent the increasing relation of 
these orbits to the horizontal one, since as we progress 
in the sequence new segments closer and closer to the 
horizontal axis are incorporated into the trajectory. 

Now, the question arises about under which conditions 
the satellite families of POs reinforce the contribution of 
the central one. As seen before, the horizontal PO satisfy 
the BS quantization condition 

TT 

kLu - -vu = 27rnH, (4) 

being nn an integer. In the same way, the BS condition 
for orbits crossing m times the surface of section reads 

fcim - —Vra = 27rn,„. (5) 

Subtracting m times Eq. to (jSJ, a new quantization 
condition emerges 

TT 

k(Lra - mLu) - -(i^m " rnvH) = 27r(n„ - mn^). (6) 

The relevant point here is that L„j — ttiLh converges 
exponentially to S^^ = —3.368 390 45 and = 
—2.991 142 22 for the two families hi and h2, respectively, 
as shown in Tabled Moreover, — rnvn is independent 
of TO, and it corresponds to t'hi = —1 and V]^^ ~ 0. In 
this way, we are in the position to assess that all satel- 
lite orbits of a given family contribute coherently when 
two BS quantization rules are simultaneously fulfilled: 
(a) the quantization of the central (horizontal) orbit [Eq. 
Q], and (b) the so-called quantization of the homoclinic 
torus according to Ref. , 

^Shi - -^i^hi = 2nn, i = l,2. (7) 



4 



Family hi 



Family I12 



3 -3.367 727 48 -2.990 915 39 

4 -3.368 367 57 -2.991 131 87 

5 -3.368 389 68 -2.991 141 81 

6 -3.368 390 43 -2.991 142 21 

7 -3.368 390 45 -2.991 142 22 



TABLE I: Exponential convergence of Lm — mLu as a func- 
tion of the number of times the POs intersect our surface of 
section for families hi and h2 of satelhte orbits corresponding 
to the primary homoclinic motion. To improve the conver- 
gence we have averaged the length of orbits m-\- and m- for 
the same family [see Fig. |21 (c)] . 
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FIG. 4: Intensity of the peak at Si^ divided by (Afc)'^ as a 
function of Afc 



It should be observed that the values of Sh^ and 
agree extremely well with those obtained numerically for 
the position of the two main peaks in the spectrum of 
Fig. 121 On the other hand, these homoclinic areas cor- 
respond in our chosen Poincare surface of section to the 
areas indicated in Fig. O (b). The arrows in Fig. O (b) 
indicate the anti-clockwise direction of the homoclinic 
motion, and for this reason homoclinic areas result neg- 
ative. 

This notorious influence of the homoclinic motion in 
the spectra, and thus in the quantum mechanics of our 
system, is remarkable, especially, when one takes into ac- 
count that the homoclinic dynamics is the origin of the 
chaotic behavior in Hamiltonian systems. In order to un- 
derstand this behavior we note that the motion has two 
components. The first one, close to the primary homo- 
clinic torus, is such that the corresponding trajectories 
interfere coherently, as it has been shown. The second 
component, associated to the non-primary homoclinic 
crosses, gives rise to more complicated self-intersecting 
tori. These excursions are much longer and the corre- 
sponding interference is much less coherent; accordingly, 
its influence in the spectra is very diluted. In this respect, 
it is more efficient to consider these long excursions, by 
the interaction with other short POs (heteroclinic con- 
nections). This point has been discussed elsewhere [Tsf . 

Finally, let us consider if the influence of the ho- 
moclinic motion on a^d survives the semiclassical limit 
?i— >0(fc = |;— >ooin our units). For this purpose we 
have performed the following calculation. Starting from 
the scaled fluctuation curve CTreh presented in the inset to 
Fig. 121 we compute the Fourier transform of the function 
for increasingly larger intervals of /cbSj {ko,ki), keeping 
fixed the value of the lower limit, fco- The results for the 
peak at are shown in Fig. 0] where it is seen that 
the corresponding Fourier intensity divided by (A/c)^ is 
an approximately constant function of the length of the 
interval, Afc = fci — fep. A similar behavior is obtained 
for the peak at 5*112. This result implies that fTioi(fcBs) 



is essentially an oscillatory function with frequencies Sh-^ 
and 5*112 ^'^'^ constant coefficients. 

In conclusion, we have gone one step further the usual 
understanding of the quantum-classical correspondence, 
by showing how homoclinic areas imprint a clear signa- 
ture in the spectra of classically chaotic systems. In par- 
ticular, we have found that the fluctuations of the relative 
widths corresponding to wavefunctions (with hyperbolic 
structure) highly localized in the vicinity of a PO oscil- 
late in a very simple way, that is clearly controlled by 
the quantization of the associated primary homoclinic 
motions. This indicates the importance of classical in- 
variants, other than those included in Gutzwiller theory 
concerning individual POs, for the systematic description 
of quantum information. 
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